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1. Overview in lieu of Introduction 



Two recent, reported results can potentially greatly enrich our understanding of quantum field 
theory (QFT). On the one hand, Komargodski and Schwimmer (KS) [I], following earlier work 
by Cappelli, D'Appollonio, Guida and Magnoli (CDGM) [2], have delineated a proof, without 
resource to perturbation theory, of an inequality satisfied when a four-dimensional QFT flows be- 
tween two fixed points of the renormalization group (RG). On the other hand, we have discovered 
closed renormalization group trajectories in four-dimensional theories, in a regime where pertur- 
bation theory is applicable pHS]. While the former result can impose restrictions on the possible 
realizations of long distance (IR) phases of QFTs, the latter gives explicitly new unexplored phases 
of QFTs. A question naturally arises as to whether these results are compatible. 

In this work we extend the work of CDGM and KS to include cycles More concretely, KS 
assume the existence of a flow from a short distance (UV) fixed point to an IR fixed point, and 
argue that the coefficient of the Euler density in the trace anomaly, denoted by a, is larger at 
the UV than the IR fixed points: aijv > °IR- This, then, is a proof of the "weak version of the 
c-theorem." We will show that we can extend the argument to include putative flows from a fixed 
point or cycle to another fixed point or cycle. Along their argument KS make several explicit and 
implicit assumptions. We make the same assumptions, and we try to be as explicit as possible 
about them. We make no additional assumptions. 

Recently, Luty, Polchinski and Rattazzi (LPR) have extended the KS argument to argue, to 
lowest order in perturbation theory, that a decreases monotonically along the RG-flow towards the 
IR [7]. This would then be a proof of the "strong version of the c-theorem." Additionally, LPR 
compute the rate of change of a as it flows along a cycle and show that a decreases monotonically 
in that case as well. Were their results correct they would render our results invalid: LPR argue 
that the monotonic decrease of a along flows on cycles rules out their existence in four-dimensional 
relativistic field theories. 

However, the computations presented by LPR are interpreted by them incorrectly. This is easy 
to see by the known fact that a increases away from trivial UV fixed points: for a Yang-Mills 
theory (YM) with beta function j3 9 = -^ g 3 /lQir 2 - fag 5 /(16tt 2 ) 2 one has [8] 

a = ao + 80Y^ + o{96) - (L1) 

Here ao is the free field theory (one- loop) value of a and ny = dim(Adj) is the number of vector 
fields. The result of LPR, that 

dt - " Cg[P j ' ° 9 ~ 32ttV' { ' 



1 More precisely "limit recursive flows." In what follows we refer to both limit cycles and limiting ergodic behavior 
simply as "cycles." 
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is incompatible with (jl.ip . We identify below the source of the error in the computation of LPR. 
We find that LPR neglect pertinent counterterms, and as a result they have in fact computed, to 
lowest order in perturbation theory, the rate of change not of a but of a different quantity, namely 
the function /3b introduced by Jack and Osborn (JO) as a candidate for a c-function. Indeed, for 
single-coupling YM models JO have already shown that /?& flows monotonically (to all orders in 
perturbation theory). 

The quantity /?& is related to a, which in JO is denoted by /3b, by 

/3 b = /3b + iw/3V. 

Here w is a function of the coupling g, and f3 g = —dg/dt is the beta function. While /3b and /3b 
agree on fixed points, the difference is parametrically large away from fixed points. In Section [2] 
we explain this in detail. There we use the method of LPR with care to include all necessary 
counterterms and rederive JO's consistency conditions, of which the monotonic flow of /3b is but 
one example. 

For models which display cycles the state of affairs is significantly more complex. In all 
these models the kinetic terms of the Lagrangian exhibit a "flavor" symmetry group Gf (that 
commutes with the gauge group). Scalar self-interactions and Yukawa couplings of scalars with 
fermions break Gf- The dependence of counterterms on the coupling constants characterizing 
these interactions is restricted by the pattern of breaking of Gf- There is a well-known, simple 
method of accounting for this. The coupling constants are treated as spurions, that is, as non- 
dynamical fields, and allowed to transform under Gf precisely so as to render the Lagrangian 
invariant under these symmetry transformations. Then, if the regulator respects the symmetry, 
so will the counterterms. It follows that the entries in the trace anomaly respect the symmetry 
too. As a is the coefficient of the G^-invariant Euler density, it is itself Gi?-invariant as well. 
And since the flow on a cycle corresponds to a Cp-transformation of the couplings, a remains 
constant on the cycle. 

This begs the following question: "how is the monotonic flow of /3b consistent with the con- 
stancy of a?" Actually, /3b is also G^-invariant, and is thus also guaranteed to be constant along 
the cyclic flow. The answer is found in the flow equation for /3b given by Osborn in [9]. His 
equation is a generalization, applicable to these more complex theories, of that found in JO. This 
flow equation is not guaranteed to give monotonic flows, but can and does give constancy of /3b 
along cycles. We review the work of Jack and Osborn concerning these more complex theories in 
Section [3j and show that f3b remains constant on cycles. We remark there that while we agree 
with JO's mathematical formulae, we disagree with their interpretation of some of their results. 
In particular, JO give a monotonic flow, but not along an RG flow, which however they interpret 
as such. 

As we have said, LPR presented an argument that gives monotonicity of the flow of a. They 
then argued that the infinite decrease of a on a cycle should be compensated by a counterterm, 
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but that in its absence it must be that no four-dimensional theory may exhibit cycles. It is a 
delicious irony that it is precisely the neglect of counterterms in the treatment of LPR that leads 
them to the wrong conclusion that the flow of a (or even /3b) is monotonic. 

In light of the lack of even a perturbative argument for the monotonic decrease of along 
RG flows, the results of CDGM and KS appear more remarkable. So in Section [5] we take on 
the task of extending the result of KS to flows between cycles or fixed points, as announced 
at the top of this introduction. We first use a non-RG flow on which a quantity Bb decreases 
monotonically, at least in perturbation theory. This quantity agrees with /3b on fixed points and 
cycles and we demonstrate that if an RG flow exists between fixed points /cycles, then the new 
non-RG flow also links those same fixed points/cycles. This non-RG flow captures the physics 
of scale-invariant theories. For example, all fixed points and cycles of the RG flow map to fixed 
points of the non-RG flow. The monotonic decrease of Bb gives still only a perturbative proof 
of a weak version of the generalized c-theorem. To go beyond perturbation theory we use the 
method of KS. The only additional assumption is that there exists a dimension-three current that 
is defined everywhere in theory space, that is, for any value of the coupling constants, such that 
it vanishes or is conserved at fixed points and generates flows on cycles. 

The natural candidate for this current was introduced formally by JO. It corresponds to a 
particular element of the Lie algebra of Gf they call S. JO showed by direct computation that 
S vanishes to three loops if the field content of the QFT consists solely of scalars, and to two 
loops if of scalars and Dirac spinors. On the other hand, the element Q of the Lie algebra of Gf 
that generates the flow along the cycle is found in our computations to arise at three-loop order 
in gauge theories that include both scalars and spinors [6]. Could it be, then, that S is non-zero 
at three loops in these theories? And if so, what is the relation between Q and 5? In Section U] 
we take on the task of computing the lowest-order contribution to JO's S for the most general 
four-dimensional QFT, compare with Q and demonstrate that S is the interpolating function we 
need to complete the extension of the KS proof of a weak version of the c-theorem to cycles. 

Let us specify here that we follow closely the notation of JO [8j, with some notable exceptions. 
From here on, following JO, in the anomaly equation we use {3 a and /3b for c and a respectively, 
although we still use the terminology c-theorem instead of the more accurate /^-theorem. Also, 
we call X a ,b,c rather than a, b, c the infinite counterterms that give rise to /3 a ,bc (having infinite 
counterterms labeled by a and c can certainly produce confusion with the corresponding "beta 
functions" that appear in the Weyl anomaly and that are commonly referred to as a and c). 
Throughout this paper RG time is defined by t = — ln(/////o), so that it increases as we flow to 
the IR. 

Many of our results are extracted from the work of JO. So it is perhaps necessary to remark 
that, besides parsing the results of JO to hopefully make them slightly more accessible to the 
general reader, we have made several novel contributions: we have discovered where in the argu- 
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merit of KS and LPR the quantity /3b is replaced by fib ( or ) m more generality, by Bb); we have 
shown that the current in JO provides a generalization of the virial away from scale-invariant 
trajectories, and is conserved at fixed points; we have extended JO's calculation of S to third 
loop order, which is the leading non- vanishing contribution, in a Yang-Mills theory with scalars 
and spinors; and finally, and most remarkably, we have used the above to prove an extension of 
the weak version of the c-theorem, that (/3(,)uv > (A>)ir, for presumed RG flows that can now 
originate or terminate on cycles as well as fixed points, valid even outside perturbation theory 
(provided the implicit assumptions in KS do not invalidate their result). 

2. Weyl consistency conditions 

In this section we review the derivation of the Weyl consistency conditions of JO. The method uses 
as a starting point the expressions of Weyl invariance used by KS and by LPR. The presentation 
is formulated so that it becomes clear that the assumptions in those works already implicitly 
lead to the JO consistency conditions. Hence, although the derivation presented here may seem 
novel, it actually follows closely JO. We have included it here for completeness, for pedagogy and 
because it makes clear that neither the results of KS nor those of LPR should be in conflict with 
those of JO. 

Let us briefly review Osborn's argument for the consistency conditions [9]. These are analogous 
to the well-known Wess-Zumino consistency conditions [10]. Let A CT VF denote the action of a 
Weyl transformation on W, the generating functional for connected renormalized Green functions. 
Because of the Abelian nature of the Weyl group, the Weyl consistency conditions follow: 

[A^A^W = 0. 

In JO the same consistency conditions are obtained by requiring finiteness of the trace of the 
stress-energy tensor in curved background and with spacetime-dependent couplings. One can also 
obtain the Weyl consistency conditions based on the arguments of LPR. 

LPR start from a quantum action So which is a function of a conformally flat metric, 7^ = 
e -2T{x) r j an( j COU pli n g constants g i (fi) (in d = 4 — e regularization, with, say, minimal subtraction 
(MS)). By rescaling the fields, which are dummy variables of integration anyway, by — > (e T ) S (f), 
where 5 is the canonical dimension of the field (as in 5 = (d — 2)/2 for scalars), and using the 
^-independence of the bare couplings, one sees that the r-dependence in So arises only due to 
the scale dependence of renormalized coupling constants, g t (e T fj,). Effectively, the regularized 
generating functional W satisfies 

W[e- 2 ^ V , u ,g\v)] = W[ V , u ,g\e^^)]. (2.1) 

Alternatively, Komargodski [TT] argues that the functional is made invariant under Weyl trans- 
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formations by adding a conformal compensator t(x). One can write 

or, equivalently, that the left-hand side is invariant under r — > r + a. Note that invariance is 
guaranteed by adding the conformal compensator and including in W all possible counterterms 
that are functions of spacetime-dependent background and coupling constants, 7^ (a?) and g l (x). 
It is from counterterms that do not vanish for spacetime-independent coupling constants that the 
P a ^-anomalies arise. It is convenient, in order to keep track of curvature-dependent terms, to 
do this in a more general background metric, 

W[e- 2 ^V(*)^(M)] = ttl7^),2V (x V)], (2.2a) 
W[e- 2 ^ lf ,Ax),g l (e-^^)} = W[ 7 ^(x), <?>)]. (2.2b) 

At the risk of restating the trivial, let us emphasize that it is not consistent to neglect the 
spacetime dependence of couplings when studying Weyl anomalies, since the Weyl transformation 
involves spacetime-dependent couplings. The counterterms associated with spacetime derivatives 
of these couplings will lead to additional anomalies. Some of these may — and as we will see, 
do — contribute to the dilaton/compensator scattering amplitude even after one takes the limit of 
flat background and spacetime-independent coupling constants. 

The approach of LPR allows one to compute quantities associated with a model in a curved 
background with spacetime-independent coupling constants in terms of corresponding quantities 
for the same model but in a flat background with, however, spacetime-dependent coupling con- 
stants. This observation is not new. For example, in JO the same observation is used precisely 
for the same purpose, namely, to compute the anomalies associated with scale transformations 
using only computations in flat space. Similarly, the approach of Komargodski allows for an 
explicit nonlinear realization of scale invariance, at the price of introducing spacetime-dependent 
coupling constants. To reiterate, in either case it is important to realize that new counterterms 
are required to render the model finite, much like counterterms involving derivatives of the metric 
need to be introduced to render finite the model in a curved background. These new counterterms 
must involve derivatives of the coupling constants and lead to new Weyl anomalies. At the end 
of this section we study how these new anomalies contribute to the Wess-Zumino action for the 
conformal compensator t(x) even after the couplings and the metric are taken to be spacetime- 
independent. For the remainder of this section we take a closer look at these counterterms, the 
anomalies they produce and the relations between them, that is, the JO consistency conditions, 
that follow from (pT2|) . 

Consider the theory in the background of an arbitrary metric Juu(x) and arbitrary spacetime- 
dependent coupling constants g l (x) corresponding to interaction terms g l (x)Oi(x) in the La- 
grangian. The arbitrary spacetime dependence of the couplings allows one to use them as sources 
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for operators in the interaction part of the Lagrangian, by taking functional derivatives of the 
generating functional with respect to g l {x). Let W stand for the renormalized generating func- 
tional. It is convenient to separate the counterterms that are independent of quantum fields from 
the rest of the action. They can be taken out of the functional integral and contribute directly 
to W: 

W = W + W c .t, 

The generating functional W results from performing the functional integral over quantum fields 
in the absence of the quantum-field-independent counterterms. The counterterms required to 
render the theory finite were first classified in JO. They consist of all possible diff-invariant 
dimension-four operators constructed out of the metric and couplings and their derivatives: 



W cX . = - J ^n~ e X-M, 

where dimensional regularization is used with d = 4 — e and 

A • & = \ a F + X b G + X C H 2 + Sid^E + l^d^d^H + \%d^d u g j G^ 

Here F is the Weyl tensor squared, G is the Euler density, H is proportional to the Ricci scalar, 
and G^ u is the Einstein tensor: 

4 2 
t-K d _ 2 K iV+ ( d _2)(d-l) K ' 

2 

^ = m _ 3)(^ _ 2) ( Rlxvp ° R V"pa ~ ^R^R^u + R ), 

1 2 

H = R, G^ v = - — - (R^v - ^lnvR)- 

The quantities above are defined as in JO, with inessential (i-dependent factors for later conve- 
nience. Each of the counterterms in A • M is an expansion in 1/e chosen to render W finite — for 
this one must in addition introduce wave- function and coupling constant counterterms, as usual. 
The coefficients A are in general functions of the couplings g l (x). 

The anomalous variation of the generating functional is dictated by these counterterms. While 
W satisfies (|2.1|) and (|2.2p . this is not true of W c .t., as can be seen by explicit computation. The 
anomaly is precisely the statement that the infinitesimal transformation r — >• r + a in (|2.2b|) . 

A ff W c .t. = Wc.t.[(l - 2<7) 7Mi „<f - oniW/dn)] - W c . t .[^ u , g% 

fails to vanish. The anomalous variation can be split into a term that would occur even if a were 
spacetime-independent plus a term proportional to the derivative of a: 

A^ anoma i y = A a W c . t . = - J ^ (*p x -M + d^a 2^). (2.4) 
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These terms again can be expanded using dimensional analysis and diff-invariance: 

Px ■ % = (3 a F + p b G + p c H 2 + xtdJd^H + ^dJ^H + \^dJd v tfG 

^ = G fiV w i d l/ g i + ^{dR) + ±i?Y^ 



(2.5) 



anc@ 



(2.6) 



+ d IM (U i V 2 g i + iV^-cU^V ) + Sijd^V 2 ^ + ^T ijk d v g i ff'^d li g k , 

up to terms with vanishing divergence. Since the anomaly gives a correction to a relation between 
finite Green functions (see 12.15]) . the anomalous terms ()2.5p and (|2,6p must themselves be finite. 

Calculations of the coefficients in j3\ ■ & and .2^ can be done using standard techniques of 
dimensional regularization with a mass-independent renormalization scheme, say MS. For now, let 
us concentrate on the relatively straightforward computation of (3\ ■ M. Since for constant a the 
transformation 5j^ u = —2a^f fiu just counts dimensions, and the dimension of the volume element 
is d while that of the operators in W c t . is four, we obtain 



(e - •& = p x -&, (2.7) 



where the beta function is 



dq i 

[i— = j3 % = -ek % g % + j3 % {g) (no sum over i). 

Here the derivative is taken holding the bare parameters fixed. k l is defined by requiring that 
the Lagrangian scales appropriately: for 0' = fj, Se (j) and g n = [i k%t g l , then «Sf ((/>', g') = /i~ e Jz? (</>, g). 
Note that f3 l di = {3 l d/dg l denotes substitution of g % by f3 % wherever g l may be found, e.g., 
^diid^g 3 ) = d^fi 3 = dif3 J d^g l , and of course respects the standard rules of differentiation. Using 
(12. 7|) . it is straightforward to show that, e.g., 



X % = (e - P k d k )s4j - MkdjP k - ^ jk dS k . 

The consistency conditions of JO can be understood as following from requiring that 
applied to the complete renormalized generating function W fails only up to the finite, anomaly 
terms. The left-hand side of (I2.2ah does not involve any counterterms from spacetime-dependent 
couplings, while the right-hand side does not involve any from a curved background. Hence, the 
counterterms in one and the other case must be related. Consider on the right-hand side of (|2.2a|) . 
for example, the counterterm 

l^iVVvV = l^ijHK^T? + • • • , (2.8) 



2 The second term involves the function of coupling constants d, which is not to be confused with d — 4 — e. Wc 
follow Osborn in this unfortunate choice of notation, hoping that with this warning no confusion will arise in what 
follows. 
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where we have expanded to lowest order in t(x). Comparing with the counterterms on the 
left-hand side of (|2.2ap . that arise solely from a curved background, we have, 

^ (X a F + X b G + X C H 2 ) = 8X b [(V 2 r) 2 - (d^rf + •••]+ 4A C [(V 2 r) 2 + •••]. (2.9) 

The X b term is a total derivative so for localized t(x) it does not contribute (recall there is an im- 
plicit spacetime integration). Matching the terms in (|2,8p and (|2.9p we find that the counterterms 
are related, 

iXc-l^pfi, (2.10) 

where the symbol ~ denotes equality up to finite terms, that is, the difference is finite as e — > 0. 
This precisely corresponds to Eq. (3.12) of JO. Applying fid/d/j, on the bare couplings to derive 
RGEs and the corresponding beta functions, one then derives from this JO's consistency condition 
(3.21a), 

8& = XijFP ~ (2.H) 

where X arises from the finite difference between the left- and right-hand sides of (|2.1U|) . and 
f3 c and xfj are beta functions for A c and ja^j, respectively. The remaining consistency conditions 
in JO can be obtained in a similar fashion. We only quote here one other consistency condition 
that plays an important role in what follows. Using (|2.2ap the lowest order terms in t{x) that 
are linear in the Einstein tensor give 

SihX,, - '//,, P. (2.12) 
If the finite difference is denoted by Wi, then one obtains 

This consistency condition is the origin of the proposal in JO for a c-function, 

h = a, + i?wi, 

which satisfies 

d^ b = l(x%^d m] )^, (2.13) 

where d^Wjj = diWj — djWi. Then its RG flow is monotonic provided the "metric" xfj is positive- 
definite, for 




To summarize, the extension ()2.2[) of the invariance requirement of LPR in (|2,ip . when applied 
to the complete set of counterterms required for finiteness when coupling constants have spacetime 
dependence, leads to the consistency conditions of JO. 
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2.1. The trace anomaly and the computation o/V^iS^ 

As formulated, the renormalized generating functional W is a function of the background met- 
ric and of renormalized spacetime-dependent coupling constants. As such we can obtain Green 
functions by functional differentiation, 



5W 1 5W 

T^(x) = -= ——- and Oi(x) = -= --7- . (2.14) 
v — 7 v \x) yj— 7 og l [x) 

These definitions are to be understood as expressions for insertions in Green functions, so that, 
for example, 0j(x) stands for the renormalized insertion of the operator which may differ from 
the operator monomial in an expectation value (JO make this distinction explicit by referring to 
the Green function insertion as [Oj(x)], but we refrain from making that distinction to streamline 
notation knowing that there is little chance of confusion). 
Using (prHD in ([2J2]) and (J23J) one obtains 

T% = ftd - fi-^ x ■ m + fi~ £ V^. (2.15) 

This is the well-known trace anomaly, accounting for the effects of curved background and 
spacetime-dependent coupling constants. However, this equation is not quite correct in the most 
generality: there are two terms missing on the right-hand side. The first is an operator that van- 
ishes by the equations of motion times the anomalous dimension of the corresponding quantum 
field. We have lost track of this term because the relation (|2.ip is only correct up to terms that 
vanish by the equations of motion. The second missing term is more subtle: we have missed 
counterterms that may be needed to render some theories finite. When the kinetic terms of 
the Lagrangian exhibit a continuous symmetry the current associated with this symmetry is a 
dimension-three operator and a new type of counterterm is required in the presence of spacetime- 
dependent couplings, that is, a counterterm proportional to the product of the current and the 
derivative of a coupling. This will be discussed extensively, and the anomaly equation will be 
fixed, in Section [3l 

Let us turn to the computation of 3?^ in (|2.4p . It follows, of course, straightforwardly from the 
definition (|2.4p . Slightly less trivial is the fact that the computation must give a finite current 
3?^. That this must be so can be seen from the trace anomaly in (|2.15p . in which all other 
terms are already finite. This means that there must be cancellations among infinite terms that 
contribute to 3fP. In fact, these cancellations are nothing but the consistency conditions, e.g., 
(|2.10p and (|2.12p . For example, the terms in fj2.4|) involving the Einstein tensor (modulo terms 
that do not vanish for spacetime-independent a) are 

J V^fG^(-8X b V II d u a-^d fl g l ^d u a) = J ^d u a 0^8^(80^ - 
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Here, in going from the first to the second line we used the finiteness condition (|2.12p and the 
definition that the finite difference is itfj. Thus we have reproduced the first term in ^ of (|2.6p . 
The remaining terms in (12. 6p can be similarly found. 



2.2. Wess-Zumino action 

The Wess-Zumino action, Wwz> is some function of t(x) that will give — AVF anoma i y upon a Weyl 
transformation, r(x) — > t{x) +a(x). Focusing on the ^-term in AWanomaiyj 

- J \f~lo~fibG-, (2.16) 

KS write the corresponding WZ term aa§| 

J ^/{r/3 b G - Af3 b [CTr,^ + r^r^vV + \(r^f] } , (2.17) 

where we have introduced the shorthand r« = d^r. However, this computation is incomplete. 
The problem with this is that we have ignored the effect of the new counterterms arising from 
spacetime dependence of the couplings. Since we will not need a Wess-Zumino action for our 
generalization of the KS argument to theories with cycles, we will not aim at being complete and 
only point out one interesting consequence here. Consider, for example, the term in 



^fd^ Wl G^d u g\ (2.18) 

Now with d^g 1 = —/3 % Tu, one has the following generalization of the Wess-Zumino dilaton action: 

J >/=7 {forG - 4(A, + \wi?) [CTr^ + t^V 2 t + £(t, m t'") 2 ] } . (2.19) 

The Weyl variation of (f2T9l) gives the sum of (|2T6|) and ([238]) (if d^g 1 = -/3V M there). The 
correction that takes into (3 b = fib + \ w ifi % is generally of lower order than the running of 
That is, ^Wi(3 l is of lower order than fib — fi^, where fi^Q stands for the free field theory value of 
fib- 

Let us be more explicit. Consider the perturbative result of JO for a YM theory with gauge 
coupling g, 

/3fc = /3b0 + 8^§F 54 + o(56) - (2 - 20) 

JO also give gw = 2ny/167r 2 + • • • , and therefore 



3 The sign in the term cubic in r is opposite to that of KS because we use the opposite sign convention for the 
conformal compensator, which gives A CT as in JO. 
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Notice however that, as already mentioned in the introduction, LPR give 

LPR fly / n n\2 

IF = ~32^ [P 1 ' 

which is incompatible with (|2.20p . However, LPR would have obtained the correct result, namely 
that they computed in fact the rate of change of fa, had they realized the dilaton scattering 
amplitude is given by the improved Wess-Zumino action in fj2. 19j) . 

It is instructive to compare the one other element of the LPR computation with the corre- 
sponding one in JO. The coefficient c g of LPR appears to play the role of the "metric" xfj m 
the consistency condition (|2,13p . In Appendix [A] we show that c g is actually related to JO's 
—2xfj- We have already encountered xtj as a metric in (|2,11|) . Explicit computations show that 
— 2x?j agrees with xfj to second order in perturbation theory in any four-dimensional theory — see 
Appendix [A] Since xtj 1S negative-definite, were xfj = — %Xij to hold to all orders in perturbation 
theory it would ensure positivity of xfj to all orders. However, this is not the case: as we show 
in Appendix [A] the equality fails, for example, at third order in a YM theory with a single gauge 
coupling. 



3. Flavor symmetries, dimension-three operators and the corrected trace anomaly 

As we have mentioned above, in deriving the trace anomaly we have missed counterterms that 
may be needed to render some theories finite. When the kinetic terms of the Lagrangian exhibit 
a continuous symmetry the current associated with this symmetry is a dimension-three operator 
and a new type of counterterm is required in the presence of spacetime-dependent couplings, that 
is, a counterterm proportional to the product of the current and the derivative of a coupling. 
Consider a theory with ns real scalar fields. The kinetic part of the Lagrangian 

exhibits a continuous symmetry under transformations of the fields 8<j) a = Wa&^b) where uj is in 
the algebra of the flavor group Gf = SO{ns). In order to renormalize this theory in the presence 
of spacetime-dependent coupling constants one must introduce a new counterterm of the form 

^c.t. = (WX^UMA, (3.1) 

with (A r j) a b = —(Nijba, that is, in the algebra of Gf- Note that this new counterterm is not 
accounted for in W c ,t. which by construction is independent of quantum fields. Note also that 
additional counterterms, symmetric under a f-> b, must also be introduced. One may integrate by 
parts to write these as terms with no derivatives acting on the quantum fields. While necessary, 
they do not play a central role in what follows. 
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To be more explicit, we consider a theory of real scalars and write for the bare Lagrangian 
= \ j^D^oaD,, 0oa + l(d- 2)4> 0a (j)o a H - ^glbcdfoafobfocfod- (3.2) 

This is written in term of bare fields (J)q. The second term is introduced to ensure conformal 
invariance of the classical action. In the potential term, the bare couplings g^ bcd are completely 
symmetric under exchange of the indices a, b, c and d. The covariant derivative, 

Dpf = (fy + A ^, 

is introduced with an eye towards including the counterterm (|3.ip . 

4)„ = A M + {D^g) I N I . (3.3) 

Here, following JO, we use the compact notation I = {abed) and we have left implicit the Lie- 
algebra indices (so that Nj = —Nj and = —A^). Note that Nj is a function of the renormal- 
ized couplings that has an e-expansion starting at order 1/e. If the theory contains gauge fields 
and some of the scalars are charged under the gauge group G g C Gf, it is straightforward to 
include an additional quantum gauge field in addition to the background field A^. 

The Lagrangian (|3,2p is explicitly locally G^-symmetric if we agree to transform the couplings 
and the gauge fields: 

Sgabcd = u ae g ebcd + permutations (5gi = (ug)i for short), 
5 An = D^u. 

The first of these is already used in defining the covariant derivative (D^g)i in (|3.3|) . It is very 
important to note at this point that if this explicit local invariance is non-anomalous it can (and 
will) be used to constrain the counterterms and the generating functional W, 

W[^(x),(n g ) I (x),QD fi n- l }=W[ 1 ^(x),gi(x),A fl ], (3.4) 

where £l(x) = exp(u(x)) G Gf- Of course, in theories without spinors the symmetry is triv- 
ially non- anomalous. Furthermore, derivatives of the generating functional with respect to the 
background field now give insertions of the scalar current. 

It is not our intention to repeat the calculations of JO in their entirety here. We will instead 
describe the main ingredients and results. We have already described the two main new ingredi- 
ents, namely, the need for new counterterms and the introduction of a background field to ensure 
invariance under Gf in (|3.4p . As before, additional quantum- field-independent counterterms are 
required. These are as in (|2.3|) but with the replacement — > to ensure Gp invariance. 
Additional counterterms involving the field strength Fn V = [Du,D u \ are also required, 

A • St = A • St + iTr(JTF^F^) + ^{^uF^D^D^j. (3.5) 
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Moreover, as advertised, new field-dependent counterterms must be included, 

=2 = rj ab (p a (j) b H + (5i) ab 4> a 4>b(D 2 g) I + |( MDtfMD^j. (3.6) 
Proceeding much as before, JO find 

T% = BjOj + ^ + {D^) T Bfa - ^f3- x ■ M + V M (^ + J» + Jg) - ((1 + 7 + S)<f>) ■ ^S . 

Many comments are in order. The last term, involving the derivative of the full action integral 
So, vanishes by the equations of motion. We have included it here for completeness. We have 
already commented that a similar term is missing from (|2,15p . The operator Oj = —j\4> a 4>b4'c4 > d 
is precisely the interaction term in the Lagrangian, and its coefficient is defined as 

^ = - {Sg) u 

where f3i = [J-dgi /dfi. The Lie-algebra element S is defined through 

S = NiBi or S = -Nj gi , (3.7) 

where Nj is the residue of the simple e-pole in Nj. Similarly, 

= p£ + D^S = pdD^g)! + D^S = P/p^)/, 

where the beta function for the background field A v is f3^ = [idA v jdyu. The beta functions for 
the field-dependent quadratic counterterms are 

The term /3t • St is the obvious generalization of (12.5P while the current is defined as in (12.6P 
but rendered G^-invariant by replacing derivatives by covariant derivatives. In addition, 3?^ has 
contributions from the new counterterms in (|3.5|) and there are additional contributions to the 
terms with and S$ that involve B^. Finally, the current is 

= (£>^) T S0, 

while the current Jq arises from the counterterms in (|3.6p and has a complicated expression in 
terms of the simple e-poles in Si and e/j (see JO for details [U Eqs. (6.21-22)]). 

For the special case of vanishing background field, spacetime-independent couplings and flat 
background the only terms that remain in the trace anomaly equation are 

2* - = B l O l - ((1 + 7 + 5)0) • ^5 . (3.8) 

We have written the divergence of the current J 11 on the left-hand side of the equation to highlight 
the similarity of this equation with that obtained when looking for scale but not conformally 
invariant theories. In that case one considers the dilatation current 

V = x v T» v - V 
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with the virial current given by 

V" = (D^fQ^, 

so that 

d^W = T% - d^V* = 0! - (Q 5 )/)C?/ - ((1 + 7 + Q)4>) ■ -^S . 

In this case Q is determined by the condition that the divergence of the dilatation current vanishes. 
This only occurs for special values of Q and the coupling constants, namely those which solve 

h - (Qg)i = o. 

It is clear that (|3.8|) suggests a relation between and V\ This relation can be easily 
understood from the behavior of QFT in curved spacetime. Indeed, it is well-known that Weyl- 
invariant theories in curved spacetime are conformal in the flat space limit [12]. The obstruction 
to lifting non-conformal scale-invariant theories to curved spacetime comes precisely from the 
existence of a non-trivial virial current. To account for the virial current, it is thus natural to 
introduce a background gauge field A^{x) which then allows for the general dilatation current 

V» = x v T» v - J*. (3.9) 

It is instructive to understand the origin of the current in JO's trace anomaly equation 
(I3.8D . The Weyl transformation for the generating functional of (12. 2p now must account for the 
transformation of the background vector field (or, alternatively, for the presence of the new 
counterterm iVj): 

W[e" Jr( * ) 7F.^W.9iW] = W[7^,A A (e T (*V), 5j (e^)/i)], (3.10a) 
W[e- 2 ^ 1 ,u,A x (e- T ^ f i),g I (e- T ^ f i)} = W[^ v , A x (jm), 9i(jj)]. (3.10b) 

Focusing on, e.g., (|3.10bp . we have indicated that all of the renormalization-point dependence 
on the left-hand side, which enters also through the counterterm in A\, carries the conformal 
compensator factor. Concentrating on the new contribution, the variation under r — > r + a gives 

A a [^(D^) T D^} = -(DrfflNihdP*)^ + ■■■ 
= -3^a J" + • • • , 

where the ellipses stand for the variations with a that are accounted for elsewhere, and we have 
used in the last step the second equation (|3.7j) . As we see, the natural outcome from a scaling 
transformation is not the canonical dilatation current x u T^ u but the current (|3.9p . 

It is now obvious, as will be explicitly shown in the next section, that corresponds to the 
virial current on cycles. More precisely, is the appropriate virial current everywhere in 
parameter space. The natural dilatation current is thus defined by (|3.9f) and differs from the 
canonical one as soon as is non-trivial. 
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The divergence of the current in (|3,8p precisely cancels the corresponding terms on the 
right-hand side, so that one has, 

=&0j-((l+ 7 )0).^Sb, (3.11) 

which is nothing but the well-known expression for the trace anomaly in flat space. This is as 
expected because (|2.2p are correct. However, in order to derive the consistency conditions in the 
most general case, one must use A3. lQf) instead of (12, 2ft to account for all possible counterterms. 
It is therefore not surprising that the work of KS led LPR to the conclusion that non-conformal 
scale-invariant theories are impossible: by neglecting these important counterterms both of them 
effectively assume that scale implies conformal invariance as a starting point. Let us turn anew 
to the JO consistency conditions. 

In the general case considered here the JO consistency conditions are modified relative to 
what has been presented in Section [2J On the one hand the conditions have to be covariant 
under transformations by the symmetry group Gf- On the other, there are additional terms that 
arise from the additional counterterms required to render the theory finite. Osborn gives the form 
of these most general consistency conditions [9|. Two conditions play a role in our discussion: 



(3.12) 



8<9//3b = xfj-Bj - B M d M wi - (d I B M )w M - (Pi9)mw m 
= x 9 u B J ~ PmQmwi - {d I f3 M )wM - {pig)uWM, 
and 

BjPi = 0. (3.13) 
In addition, covariance under Gf gives, e.g., 

(ug)idiP b = and (ug)idjS = [u, S]. (3.14) 

Of course, the first of these applies to any Gj?-invariant while the second to any antisymmetric 
tensor (for example, any Lie-algebra valued function). Using the first of (|3.14p in (|3.12p gives a 
nontrivial relation among several beta functions: 

{ug)i [XfjBj - B M d M wi - {diB M )w M - (Pi9)mw m ] = 0, (3.15) 

or, equivalently, 

(ug)i [x 9 u B J ~ PmQmwi - (9iPm)w m - {pi9)mwm} = 0. (3.16) 

These conditions can be used to understand aspects of the flow of Consider the flow 
defined by some arbitrary function fi(g), 

f - -/'«-»• 
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If one takes fj = fij then the flow can be identified with the RG flow, with rj = t = — ln(/////o). 
From (|3.12p we have 

where Bb = fib + ^Bjwj, and 

dfib 

~ 8 ~drj~ = X ?^ /jBj ~ flPMd[ M wi\ ~ {pig)MfiWM- 

Three special cases are of most interest. Consider first //((?) = (ojg)j. From the second equation 
in (|3.14|) we see that on the flow S is constant. This is a recursive flow (cycle or ergodic). It 
follows from the G,F-invariance of Bb and fib that these remain constant on the flow. This is a 
consequence of the detailed cancellations that must be satisfied by the beta functions in (|3.15|) 
and (|3.16|) . This general result can be applied to actual limit cycles, for which u = Q. Based on 
an incorrect equation, roughly (|3. 12j) with 5, wi and pj set to zero, LPR argue that cycles cannot 
occur in four-dimensional theories. They argue that fib decreases monotonically, dfib/dt < on 
the cycle, which is inconsistent since there is no counterterm to remove the infinity as fib flows 
to the IR. We now see that indeed there are counterterms, missed by LPR but known to JO, 
that not only account for the would-be divergence but in fact ensure Gi?-covariance and hence 
constancy of fib (and fib) on the cycle. 

The second special flow we consider is that of // = fij, that is, actual RG-flow. Now we have 

- 8^ = X 9 ijfiiB.l ~ {pig)Mfiiw M - (3.17) 

If we were to ignore the last term on the right-hand side we could establish a perturbative c- 
theorem for fib- Indeed, to two loops Bj = fij and Xij = ~^X a ij > so the right-hand side of ()3.17|) 
would be positive-definite along a perturbative flow. However, the last term is parametrically of 
the same order as the first on the right-hand side of (|3,17p precluding a perturbative proof of the 
strong version of the c-theorem. 

The third and last special case we wish to consider is that of // = Bj. In that case, 

-8^ = xf J 5/ J Bj, (3-18) 

where we have used the consistency condition (|3. 13j) . This shows that there are flows on which 
some quantities decrease monotonically. However, these flows are not RG flows, even though 
they are perturbatively very close to RG flows. If one could establish that these flows have as 
endpoints the same fixed points (or limit cycles) as RG flows one would have a proof of the weak 
version of the c-theorem. We return to this question in Section [5j 

It is worth noting that J^, although not conserved, has no anomalous dimension and thus 
dimension exactly three on a cycle. This observation is important because it contradicts the view, 
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espoused by Osborn, that one can take Bj to be the beta function. If this were the case one 
would have conformal invariance on cycles (which would then be interpreted as fixed points), 
which however is impossible in light of unitarity and the existence of the non-vanishing non- 
conserved current 3^ of dimension three. It is clear that the relation between /3/ and Bj is not 
simply a change of renormalization scheme, since some fixed points of Bj are not fixed points 
of For this reason, in the rest of this work we will continue to make explicit the distinction 
between RG-flows and tj-flows generated by Bj. 

4. S is Q (on cycles) 

We alluded above to the similar role played by the current J M that arises in the trace anomaly 
and the virial current that defines a conserved dilatation current on a cycle. In this section we 
investigate the relation between them, or more specifically, the relation between the Lie-algebra 
valued function S and the Lie-algebra constant Q defined only on the limit cycle. 

What is known about S? JO have shown, by direct calculation, and we have verified, that 
in a scalar field theory S vanishes up to third order in the loop expansion. The result holds 
even if gauge fields are included and the scalars are charged under the gauge group. For theories 
with scalars and fermions, JO have shown, and we have verified, that S remains zero to two 
loops. However, this is consistent with a possible equality of S and Q on cycles. Indeed, we have 
obtained previously that Q is of third order in the loop expansion in Yang-Mills theories with 
scalars and fermions, while in purely scalar theories we were unable to find a non- vanishing Q up 
to third order in the loop expansion. 

As expected from the discussion above, we will show that (up to conserved current) 

1. S is Q on cycles, 

2. S vanishes at fixed points. 

In light of these results the computation of Q can be tremendously simplified given an explicit 
expression for S. The procedure to determine Q involves determining first the beta functions for 
the coupling constants to second order in the loop expansion for scalar self-couplings, to third 
order in the loop expansion for Yukawa couplings and to fourth order in the loop expansion for 
Yang-Mills couplings, and then solving the system of nonlinear coupled equations f3 g = and 
Pi = {Qg)i ( we implicitly use here that gj can also stand for Yukawa couplings). Since S must 
have a perturbative expansion that starts at third order in the loop expansion, to determine Q 
from S it suffices to evaluate it with coupling constants on the cycle computed to lowest order 
in the loop expansion. So Q is obtained from S by determining the zeroes of the one-loop beta 
functions (two- loop for gauge couplings): if S = on the zero of the beta functions, the zero is a 
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fixed point of the RGE, but if S 7^ on the zero, then the zero is a point on a cycle and Q = S 
there. 

To this end an explicit, three-loop expression for S is required. But as pointed out above, 
there has been no computation of S to the order where one would expect it to be non-vanishing 
if S were to equal Q on cycles. We have endeavored to compute S to third order in the loop 
expansion for a general theory containing ns real scalars and n/ Weyl spinors, possibly charged 
under a gauge group. The potential in the Lagrangian is 

-V = %Kbcd</>a</>b<f>c<l>d + (^y a \ij<Pa^i1pj + h.C.). 

The details of the computation are spelled out in Appendix [Bj The surprisingly simple result is 

(16vr 2 ) 3 5 a6 = § ti (y a y*y d y*)X bcde + § tv{y a y* c y d y d y b y* c ) - {a <-> b} + h.c. 

We have evaluated this expression on the fixed points and cycles of the theories we explored 
in [31 [U [6] and found that in each case S vanishes at all fixed points and equals our previous 
determination of Q on all cycles. It should be mentioned that the models in [3l[5] are in d = 4 — e 
dimensions, with e small but positive. In these cases the consistency equations still apply provided 
one carefully tracks the presence of e in the beta functions, particularly in j3 % = —ek % g % + f3 l and 
correspondingly in B l . 

Now for the proof of the propositions above. First we show that S = Q on cycles. Consider 
the rj flow with fj = Bj, with boundary condition that at rj = the point gi(0) is on the cycle. 
Then -6/(0) = ([Q — S]g(0))i is in the Lie algebra of Gp and the left-hand side of (|3.18j) vanishes. 
Since Xjj = ~^Xu to second order in the loop expansion, and — is positive-definite, (|3.18p 
gives 5/(0) = 0. This implies S = Q + AQ on cycles, where (AQ5)/ = 0. But if AQ ^ this 
corresponds to a conserved current, and we are free to redefine the scale current by a conserved 
current, Q — > Q + AQ. Hence, S = Q on cycles. 

For theories with two scalars there is an alternative, perhaps simpler proof that S equals Q 
when evaluated on a cycle. Consider (|3.18|) evaluated on a point on the cycle. It is easy to 
show that 5 is a constant on the cycle: —dS/dt = fiidiS = (Qg)idiS = [Q,S] = 0, where the 
last two steps follow from f)3. 14|) and the fact that, for two flavors, the flavor group, SO(2), 
is Abelian. Now, as before, we consider the r/-flow defined by the Bj function starting from a 
point on the RG-cycle (we make the distinction of the actual RG-cycle and a ry-cycle explicit, 
to avoid confusion). The flow is defined by —dg/drj = Bj = f3j — (Sg)j = ([Q — S]g)i, where 
the last step follows from assuming the initial point is on the RG-cycle and then noting that the 
solution corresponds to a trajectory that traverses the same cycle but at a different angular speed 
(the angular speeds are Q12 and Qu — S12 for the RG- and recycles, respectively). Therefore 
the 77-cycle is generated by a trajectory in Gf and it follows that, just as for an RG-cycle, any 
Gi?-invariant remains constant on the ry-cycle. But the consistency condition (|3.18p then implies 
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that Bj = ([Q — S]g)j = on the cycle. Since Q and S are each characterized by a single number 
the only solution is S12 = Q12 (on the cycle). 

It is easy to show that (Sg)i = at a fixed point, and this is consistent with the notion 
that a fixed point corresponds to the case (Qg)i = 0. To see this, notice that at a fixed point 
Bj = —(Sg)j so at that point the flow corresponds to a first order ^^-transformation. That is, 
the first derivative with respect to 77 of GF-invariants vanishes at the fixed point. Hence, (13.18|) 
gives that x 9 ij{Sg)i{Sg),j = and hence (Sg)i = at the fixed point. The solution is that either 
S = at the fixed point, or there is an emergent symmetry at the fixed point, and is the 
corresponding conserved current. This completes the proof of the two propositions above. 



5. An improved c-theorem 

As we have seen, the consistency relations of JO lead to a monotonic tj-fiow in perturbation 
theory, 

= BjdjB b = Ix'jjBiBj > 0, 

where we have defined Bf, = (3^ + ^Bjwj. Only the last step in this sequence of relations invokes 
perturbation theory, for the positivity of the metric %% is established perturbatively; see Appendix 
lAl However, generally this 77-flow does not correspond to an RG-flow, so one cannot conclude from 
this that the strong form of the c-theorem — that c = Bb decreases monotonically along the RG- 
flow — applies. One could try to use instead the consistency condition for the RG-flow of /3f, in 
(|3.17p . but the right-hand side of that condition is not necessarily positive-definite. Alternatively 
we use the fact that Bj = at fixed points and on cycles, and show that there are 77-flows 
between fixed points and cycles (or between fixed points and fixed points or between cycles and 
cycles) and hence prove the weak version of the c-theorem, that is, that (/3fe)uv > (A))ir, where 
the subscript refers to fi b evaluated at the corresponding end-point (or end-limit-cycle) of the RG 
trajectory. 

We assume there is an RG flow with endpoints on fixed points or limit cycles, 

with gi(— 00) = (gi)\j\r and 5/(00) = (gi)iR. In the case of limit cycles the end points refer rather 
to limit trajectories. Now, consider the solution to the 77 flow, 




= BMv)). 



It is given by 

g(rj) = F(rj)g(ri), 
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where 

F(rj) = T ( exp - f dr/S(r/) 



OG 



is an element of Gf and T is the 77-ordered product. As such, /?/(<?) = Pi{Fg) = (Ff3)j(g) so that 
the fixed points and cycles (gi)uv and (gi)m are mapped to physically equivalent fixed points 
and cycles (gi)xjv and {gi)m- Since is (^-invariant and Bj = on fixed points and cycles, it 
follows from the monotonic flow of Bf, that ((3b)\JV > (A>)lR- We emphasize that this proof relies 
on perturbation theory to establish the positivity of the metric xf j- For a non-perturbative proof 
we turn to the method of KS. 

Let us review the argument of KS. Our presentation is closer in spirit to that of LPR. We 
will try to note explicitly when implicit assumptions in that argument are made. While plausible, 
these assumptions should be justified for the theorem to be established. We deviate from both 
presentations in that we do not derive nor use a Wess-Zumino dilaton action for, as we will see, 
this is not necessary for the computation. 

Consider the four point function of the operator ^d^T?^ = ^d^{x v T^ v — J M ) in an arbitrary four- 
dimensional theory which is classically scale- invariant. In the absence of dimension-three currents 
this is just the four point function of ^T^, but we can more generally consider in one swoop 
the case where dimension-three currents are present. Furthermore, we will consider kinematics 
such that pf = 0, i = 1, . . . ,4, for the momenta pi of the four insertions. This is equivalent to 
computing the t(x) scattering amplitude A(s, t) with the on-shell condition V 2 r = 0, so that the 
Mandelstam variables satisfy s + t + u = 0. 

Now, we will assume that the forward scattering amplitude «4f w( j(s) = {s,0) exists, that is that 
the limit t — > of (s,t) exists. This could fail if A(s,t) had terms of the form, e.g., ~ s 2 ln(t). 
As we have explained, the response of the theory to Weyl rescaling is measured in the general 
case by d^V 1 , so A{ w d{s) can be computed by taking four r(x)-derivatives of the generating 
functional and then taking the metric as flat, the coupling constants to be spacetime-independent 
and the background field and the conformal compensator to vanish. Alternatively, and more 
straightforwardly, one can work with a conformally flat metric and having the only spacetime 
dependence in couplings and A^ arise through the dependence on the conformal factor t(x), so 
that one merely needs to take r(x) = after four times differentiating W. Now the first derivative 
simply gives the dilatation anomaly equation 

= B I 1 + + (D^fBfo - ^(3- x ■ % + + Jg) - ((1 + 7 + S)<P) ■ ^S . 

One need only thrice differentiate this equation to obtain the four point amplitude of V M X> M . Note 
that on fixed points and cycles, where we will need this, the first term vanishes since Bj = 0. Also, 
the last term, which vanishes by the equations of motion, can be ignored for the computation of 
the amplitude. Most of the remaining terms vanish once the couplings are taken to be spacetime- 
independent (and the metric flat and A^ = 0). The remaining terms arise from the G and H 2 
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terms in /3\ ■ ffl. For a conformally flat metric, 7^ = exp(— 2t(x))t]^ u , one has (in d spacetime 
dimensions) 

e -^G = 8(d 2 r) 2 - 8 W >^ - 16t, mV t^ - 8(d - 3)t^8 2 t + 2(d - l)(d - 4)( V r^) 2 

The cubic term in H 2 vanishes for an "on-shell" conformal factor d 2 T = and so the only 
contribution to the "on-shell" forward scattering amplitude is from G: 

«4fwd(s)| F P or cycle = -pb{s 2 + t 2 + U 2 )\ t=Q = -2(3 b S 2 . (5.1) 

We assume there exists an RG trajectory from a UV fixed point or cycle to an IR fixed point 
or cycle. On this trajectory this equation no longer holds. However, we can inspect limiting 
behavior. Since -4f w d/s 2 depends on s only through the dimensionless ratio /i 2 /s, its behavior is 
dictated by the renormalization group. Hence, 



lim 



lim 

s— ¥00 



^fwd ( S ) | FP or cycle 



-2(A,)uv 



and 



,. -4fwd(s) «4fwd( s )|FP or cycle n/oN 

hm — = hm 3 - — = -2(/3 b ) m , 



s->0 



where (/3b)uv an d (/3&)m are the limiting UV and IR values of /3 b along the trajectory and 
correspond to those on the fixed point or cycle. LPR do a careful analysis to validate the 
approach to these limiting values using conformal perturbation theory. 

Following LPR we next consider the integral of A{ w( i(s)/s 3 over the contour in Fig. [TJ The 




Fig. 1: The contour of integration for 



ds 



-4fwd(s)- 



integral over the semicircle I\ cannot be easily computed, but in the limit that the radius of the 
semicircle vanishes it is reasonable that one can use the limiting value, 



ds , . 
- o'>Afwdv*J 



f dn 

/ -2(/3 b ) m = 2m((3 b ) m , 
Jh s 



(5.2) 
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where the last step corresponds to taking the vanishing limit of the radius of the semicircle I\. 
Similarly, the large circle I3 gives 

r ds f ds 

/ -T^fwdW« / — 2(/3 6 ) UV = -27rz(/3 b ) uv . (5.3) 
J Is s Jh s 

It follows from Cauchy's theorem that 

1 f ds 

(A>)lJV - (A>)lR = J i ^3^fwd(s) 

ds 

-^lm(A iwd (s + i0)), 

where in the last line LPR assume crossing symmetry to write -4f w d(— s + iO) = Al- wd (s + iO). 
Finally, the KS argument invokes the optical theorem that relates the imaginary part of the 
forward scattering amplitude to a positive-definite cross section to conclude that 

(AOuv - (ft)m > 0. 

We note in passing that the optical theorem is known to apply for forward scattering amplitudes 
of (on-shell) physical particles. It is not clear a priori that it applies to Green functions of 
composite operators at pf = 0, even if it corresponds to the scattering amplitude of would-be 
dilaton scattering. We think the assumption of positivity is reasonable, so we press on. 

What steps in the argument above require special attention when the theory admits dimension- 
three currents? As we have pointed out, the dilatation current now has an additional term, 
but we have already accounted for this in the presentation above. Throughout the flow this 
makes no difference to the argument above, since the positivity of the integral over the segments 
^3 of the contour follows from the optical theorem. For cycles one is not free to ignore the t{x) 
dependence of the couplings or the background vector field in the anomaly equation. But on 
the cycle the couplings are covariantly constant. Hence, the terms that vanish at fixed points 
because of the constancy of couplings also vanish for cycles, but now because they are covariantly 
constant. Finally, the validity of the limits in (|5.2p and (|5.3p need to be established anew for 
limit cycles. However, the same method of conformal perturbation theory may now be applied, 
using instead scale-but-not-conformal perturbation theory to establish the result. Since it is only 
scaling that is used in this step of the argument by LPR, the proof goes through as presented 
there. 

6. Conclusions 

We have shown that a weak version of the c-theorem can be shown along the lines of the arguments 
presented by Cappelli, D'Appollonio, Guida and Magnoli and by Komargodski and Schwimmer 
for the more general case that a renormalization group flow goes from a fixed point or cycle to 
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another fixed point or cycle. To this end we demonstrated that the quantity that plays the role 
of c is Jack and Osborn's B^ which is closely related to the a-anomaly (fib in the notation of Jack 
and Osborn); these quantities agree at fixed points, but are not generally the same. We showed 
that given an RG trajectory linking fixed points and/or cycles, there is a non-RG trajectory 
linking them too on which Bf, decreases monotonically from the UV to the IR. 

We showed where the arguments presented by Luty, Polchinski and Rattazzi, that no renorma- 
lization-group limit cycles can exist in four-dimensional theories, are in error. In contrast with 
their results, we showed that /3& is constant on a cycle. To this end we reviewed (and presented 
a somewhat different derivation of) the consistency conditions of Jack and Osborn. The crucial 
observation is that Luty et al neglect to include the effect of essential counterterms. 

We presented a calculation of the Lie-algebra function of coupling constants S introduced by 
Jack and Osborn. This is the first calculation of S to an order (third) in the loop expansion where 
it does not vanish. We then proved that S = on fixed points and that S precisely corresponds 
to the generator Q of limit cycles when evaluated at any point on the limit cycle. This is a major 
improvement on the method of searching for limit cycles: one merely needs to find zeroes of the 
beta functions to the first order in the loop expansion (second order for Yang-Mills couplings) 
and evaluate S there. If S = the zero corresponds to a fixed point, while if S = Q 7^ the zero 
corresponds to a limit cycle with Q the generator of the cycle. 

Some questions remain which we intend to turn to in the future. Among them are: 

• Are there renormalization group flows between fixed points and cycles? 

• Are there models in four dimensions with limit cycles having tree-level bounded scalar 
potential? 

• Are there limit cycles in super symmetric four-dimensional field theories? 
We look forward to answering these questions. 
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Appendix A. The relation between the metrics xfj anf i X% 

The coefficient c g of LPR appears to play the role of the "metric" xfj i n the consistency condition 
(|2.13|) . As we mention in the end of Section 12.21 and elaborate on further here, this is not the 
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case. To see the connection with the work of JO, following LPR we write 

A^ anoma i y = ^ J d i xd i yT(x)T(y)(@(x)e(y)}, 
where = (3 l Oi, and therefore 

^AVK anomaly = I J d 4 xd 4 yT(x)T(y)j t (Q(x)Q(y)). (A.l) 
In Ref. [9l Eq. (3.18b)] Osborn finds the RGE for the product of two local renormalized operators, 

The quantity xtj can be thought of as the beta function associated with the counterterm needed 
in order to renormalize the correlator (Oi(x)Oj(O)). Because of this, it is negative-definite. Now 
since —d/3 l /dt = I3 3 djj3' l 1 it is easy to see that 

|<©(x)©(0)> = x%W j d 2 dH^{x). (A.2) 

Using this in (|A. 1[) we see that the metric of LPR is —2x a , which is always positive. This suggests 
the question: is there a relation between x 9 an d X a ? 

In the specific example of a gauge theory with a simple gauge group G and charged fermions 
in some representation, JO give Eqs. (5.12)], at two loops, 

-,2 



a 1 q n V 
v = v y = 

X 2 X 8vrV 



I + i 17C - y R I h 



^ifc' ( A - 3 ) 



where t r (£adj^adj) = C5 ah , R 18 similarly defined for the representation of the fermions, and ny = 
dim(Adj) is the number of vectors. However, the relation x 9 = ~2x a °f (|A.3|) does not hold in 
general, and so the task of computing x 9 is complicated. Nevertheless, Weyl consistency conditions 
give the general relation between x a an d X 9 13 Eq. (3.23)]: 

Xfj + IXij ~ Xi 3 kP k = ~P k d k V l0 - drfi V kj - d 3 (3 k V ik , (A.4) 

where 

^ij = Xijk' 



C y ij = X a ak kk 9 k (no sum over k), 



and 

(.ijk u kXij ~ 2~KXikj 



Xijk &kXij 2^-Xikj Xjki) ; 



with Xijk necessary to regulate infinities in three-point functions, and ( defined as an operator 
counting the number of loops, whose form can be read off from 0(e) terms of the finiteness 
condition (3.9e) of JO: 

C,Vij = (1 + k k g k dk)Vij + 2k l Vij (no sum over the index carried by k), 
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(cf. [8, Eq. (3.16b)]). 

In our gauge-theory example we have 



dV 80 9 
X 9 + 2 X a - x a 9 = -P'jt ~ 

og dg 



(V = \x a g, 



(A.5) 



where (V = (2 + ^gd/dg)V = (2 + hd/dh)V, the beta function for the gauge coupling is 



V 



-Poh- 0ih 2 + 0(h 3 ), 0o = -(11C7 -312), 



and 



where x b is given at two loops by 



X 



6(2) 



X" 



ny 



d X a b 



It follows that 



Expanding V = Vq + v\h + 



X 



4.V i + - 3 (nc-4R)h 



^—[(llC-4R)h + 0(h 2 )}. 



37r 2 g 3 

gives QV = 2fo + 3v\h + 
ny 



\gx a , or 



With this, gives 



V 



X 9 + 2x a 



1927T 4 

n v o 



96vr 4 



(11C-4R)+0(/j). 
(llC-412)/i + 0(^ 2 ) 



(A.6) 



and, therefore, beyond two-loop order, \ g 7^ ~^X a - 

To summarize, as shown in JO the negative-definiteness of xtj follows from the two point 
function of the operators Oi and holds generally. In general xfj + ^Xfj f aus to vanish beyond 
the first few orders in the loop expansion. The positivity of x% may fail non-perturbatively (for 
example, if its perturbative expression has finite radius of convergence). 



Appendix B. How to calculate Nj and S 

The calculation of JO's Nj proceeds order by order in perturbation theory. In this appendix we 
calculate contributions to Nj in a quantum field theory with real scalars and Weyl spinors up to 
two loops, and we also perform a three-loop calculation of the part of Nj that is needed in order 
to compute S. 

As can be seen from (|3.ip . in order to calculate Nj we need to compute self-energies of scalars 
but with coupling constants as spectator fields. Equivalently, the calculation can be done by 
considering scalar self-energy diagrams and letting momentum come in from external legs and go 
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out through couplings. From these diagrams we can then pick up the contribution linear in the 
momentum of the field and linear in the momentum of the coupling. After we antisymmetrize, 
we have a contribution to Nj. 

It is perhaps helpful to remind the reader here that in a theory with scalars and fermions the 
/ index can be either (abed) or (a\ij). Let us also remark that S appears first at three loops in a 
theory with scalars and spinors. The reason is easily seen from (13. 7|) : a diagram that contributes 
to N will only contribute to S if it is not symmetric under a •<->■ b. As it turns out there are no 
such diagrams in scalar self-energies at one and two loops, but there are two such diagrams at 
three loops. Consequently, even if the theory contains gauge fields, diagrams with gauge fields 
will not contribute to S at three loops, but certainly will do so at higher order. Therefore, even 
in a gauge theory we don't need to include gauge fields in our leading-order calculation of S. 

B.l. One loop 

At one loop the calculation proceeds with no subtleties since renormalization is straightforward. 
The two diagrams that contribute to Nj and their corresponding counterterms are shown in Fig. [2J 




Fig. 2: Diagrams that contribute to N a uj at one loop and their corresponding counterterms. 

A straightforward calculation gives 

(N clij ) ab = -jjr^ \{v*a\i^bc ~ y* m o~ac) + finite, 

and there is of course a complex conjugate (iVV.)^. 

In order to simplify the notation we write the result for the residue of the simple e-pole in Ni 
in the form 

16ir 2 (Nj) ab d» gi = -i[tv(y a d"y* b ) + h.c. - {a b}], 

where gi on the left-hand side stands here for y c ^j or y*uj- Selecting the appropriate derivatives 
one easily reads off the corresponding Nj. Our result reproduces JO's equation (7.16) for pi when 
we use Dirac spinors. 

B.2. Two loops 

At two loops there are three Feynman diagrams that contribute to Nj, listed in Fig. [3J The 
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Fig. 3: Feynman diagrams that contribute to Nj at two loops. 



calculation of the residues of the simple e-poles of Nj requires now a subtraction of sub divergences, 
something that proceeds, for the most part, in the usual way. However, there is a small subtlety, 
not seen in the usual treatments of renormalization, that we would like to point out. Clearly, the 
two right-most diagrams of Fig. [3] have sub divergences so we have to add to them the diagrams 
with the insertions of the corresponding counterterms. For the right-most diagram the graph with 
the insertion of the counterterm is 




Now, when the momentum that comes in from, say the left external leg, flows out through the 
counterterm, then there are two diagrams that contribute, namely 



P 




and 




where the momentum exits to the north-east or to the north-west depending on which vertex it 
flows out of in the original diagram in Fig. [3l In both cases the counterterm is the same, but 
the diagram with the insertion of the counterterm is different as a result of the difference in the 
momentum of the internal leg that the counterterm picks up. That is, had we retained different 
momenta for the various vertices, there would be two momenta associated with the counterterm. 
The two loop result for Nj, previously unpublished, is 

(Wn 2 ) 2 (N}) ab d" gi = - ±\ acde d»\ bcde + [\tx{y* a y c d^yty c ) + § tr(y* a d»y c y* c y b ) 
+ i HylVc^ytyb) + I tr(y*y c y*^2/ 6 ) + h.c] - {a b}. 

It follows that S vanishes at this order. This can be seen, term by term (when anti-symmetrized 
in a and b) by replacing gi for d^gj. 
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B.3. Three loops 

At three loops there are many diagrams that contribute to Nj, but only four are not symmetric 
under a f-> b and thus end up contributing to S. These diagrams are shown in Fig. [H and we 
here only compute their contributions to Nj. 




Fig. 4: Diagrams that contribute to Nj not symmetric under a -H- b and hence contribute to S 
at three loops. 

From these diagrams (and the corresponding counterterms) , using the methods for the calculation 
of three-loop integrals given in [13], we find 

(levr 2 ) 3 ^ 1 ),^/ D - ±tT(y a d»y* c y d y* e )\ bcde - ± tT(y a y*d^y d y* e )X bcde - \ tv(y a y* c y d d^y* e )X bcde 

~ ^^(yaVcydytW hcde - ^tr{y b d^y*y d y* e )X acde - ^ tr(y 6 y*d M y d y*)A aC(fe 

- ^ti(y b y*y d d^y* e )Xacde - ^ tr^ y b y* c y d y* e )X acde - ^ tr(y a d^y*y d y d y b y* c ) 

- h iT ^yay*c d ^ydydVby*c) - §^^2/*^^^*) - h iT ^ay*cydyWy b y* c ) 

- ^tr(y a 2/*y d y^ 6 ^y*) + ^ti{y a d^y* c y c y d y b y d ) - ^tr{y a y* c d^y c y* d y b y d ) 
~ ®tr(y y*y c 5 /i y^ 6 y^) - ^ ti{y a y* c y c y d d^y b y* d ) + i tr(y a y*y c y^y 6 d^) 
+ h.c. - {a 6}, 

and since 

S 1 = -kiNjgt = -Nl bcd X abcd - {jN^yaiij + h.c.) 

we finally obtain 

(16vr 2 ) 3 5 afe = ^tr(y a y*y d y*)X bcde + § tr (y a y* c y d y* d y b y* c ) + h.c. - {a O 6}. 

As already remarked in the main body, evaluating this on points in coupling space where we have 
found fixed points and cycles in Refs. [SIEIE], we find that S vanishes at all fixed points and 
equals Q on all cycles. 
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